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ABSTRACT 



[T^ , Context. The COSMOS survey is a multiwavelength survey aimed to study the evolution of galaxies, AGN and large scale structures. Within this 
HH ■ survey XMM-COSMOS a powerful tool to detect AGN and galaxy clusters. The XMM-COSMOS is a deep X-ray survey over the full 2 deg 2 of 
^"J the COSMOS area. It consists of 55 XMM-Newton pointings for a total exposure of~1.5 Ms with an average vignetting-corrected depth of 40 ks 
, across the field of view and a sky coverage of 2. 13 deg 2 . 
O ,■ Aims. We present the catalogue of point-like X-ray sources detected with the EPIC CCD cameras, the logN-logS relations and the X-ray colour- 
I • colour diagrams. 

O ' Methods. The analysis was performed using the XMM-SAS data analysis package in the 0.5-2 keV, 2-10 keV and 5-10 keV energy bands. 

| Source detection has been performed using a maximum likelihood technique especially designed for raster scan surveys. The completeness of the 
CZ2 , catalogue as well as logN-logS and source density maps have been calibrated using Monte Carlo simulations. 

Results. The catalogs contains a total of 1 887 unique sources detected in at least one band with likelihood parameter det_ml> 10. The survey, which 
shows unprecedented homogeneity, has a flux limit of ~1.7xl0~' 5 erg cirT 2 s , ~9.3xlO~ 15 erg cirT 2 s -1 and ~1.3xl0 -14 erg cirT 2 s _1 over 90% 
of the area (1.92 deg 2 ) in the 0.5-2 keV, 2-10 keV and 5-10 keV energy band, respectively. Thanks to the rather homogeneous exposure over a 
^ , large area, the derived logN-logS relations are very well determined over the flux range sampled by XMM-COSMOS. These relations have been 



X 



compared with XRB synthesis models, which reproduce the observations with an agreement of ~10% in the 5-10 keV and 2-10 keV band, while 
in the 0.5-2 keV band the agreement is of the order of ~20%. The hard X-ray colors confirmed that the majority of the extragalactic sources in a 
| bright subsample are actually Type I or Type II AGN. About 20% of the sources have a X-ray luminosity typical of AGN (L x > 10 42 erg/s) although 
, they do not show any clear signature of nuclear activity in the optical spectrum. 

I , Key words. Galaxies: active, Cosmology:large-scale structure of Universe, X-rays: diffuse background, X-rays: galaxies 

o; 

o\ . 

© ■ 1. Introduction lElvis et all d2009l) (hereinafter C-COSMOS). 

_ _ . , . ^ In this paper we present the X-ray pointlike source catalogue of 

The Cosmic evolution survey (COSMOS, | Scoville et al. || 2007| ) ^ 1>5 Mfj XM M-COSMOS survey together with the observa- 

with its 2 deg of multiwavelength data is an exceptional lab- tion ^ data ductSi i ogN _i ogS relations and colour-colour 
; oratory to study Active Galactic Nuclei (AGN) galaxies large lots A sub le of the first year of XMM-COSMOS data has 

scale structures of the Universe and their co-evolution. The sur- beefl ented in C a PP elluti et al. (2007a) (hereafter Paper II) 

vey uses multi-wavelength imaging and spectroscopy from X- togeth er with a detailed overview ofthe data analysis techniques, 

ray to radio wavelengths, including HST, Spitzer and GALEX Herg wg t data of ^ ±e observi les> with im p rove d 

imaging. The size ofthe survey has been chosen to sample large- source itioni hi her counti statistics and more ise 

scale structures with linear sizes of ~50 Mpc h at z=l with ^ nhotometrv 

highly reduced 'cosmic' or sample variance. 0ptical identifications of XMM-COSMOS sources, performed 

During the A03-A04 and A06 cycles XMM-Newton (surveyed b ^ advantage f the precise source positioning achieved 

2.13 deg* of sky in the COSMOS field m the 0.5-10 keV en w ^ ththe L mple J ntarvCha V a observations, will be presented 

ergy band. The total exposure was -1.5 Ms split over 55 EPIC fa anQther (Bmsa et al 2QQ9) 

pointings. The average resulting exposure over the field of view The combination of the moderate i y deep flux limit and the wide 

is -68 ks. The central 0.9 deg* of the COSMOS field also has effectiye ^ (flux Hmit of ^ 1>7xl0 -is cm -2 g -i in the q.5-2 

been observed in X-rays w ith Chandra for a total of 1 .8 Ms by kgV band oyer { m deg2) of the XMM-COSMOS made possi- 

c , a- . . . . XT ^ „ . ■ ble the compilation of a sample of sources with low influence of 

Send offprint requests to: N. Cappelluti , \ , r . , . . , . . 

* T3 u »■ i . ■ i ..i VA ,» . \j . cca- the so called sample or cosmic variance. Indeed, in paper II, 

Based on observations obtained with XMM-Newton, an ESA sci- r • 

ence mission with instruments and contributions directly funded by assuming these survey parameters, we estimated that in XMM- 

ESA Member States and NASA; also based on data collected at the COSMOS the fluctuations of the source density due to cosmic 

Canada-France-Hawaii Telescope operated by the National Research variance are <5%. Furthermore, the tiling of the observations 

Council of Canada, the Centre National de la Recherche Scientifique was chosen to maximize the uniformity of the sensitivity over a 
de France and the University of Hawaii. 
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large area of the field. 

These particular characteristics, together with the multitude of 
multiwavelength information available, were designed ad hoc to 
study the large scale structures traced by X-ray emitting objects 
like AGN and galaxy clusters and their co-evolution (see e.g. 
Cappi et a ll200lHCappelluti et al.l2005l,l2007bHBranchesi et all 
2007; iKocevski et al.l 120081) . In addition these characteristics 
make the survey sensitive enough to study the evolution of super- 
massive black holes in the Universe up to high-z. Considering 
the high throughput of XMM-Newton at high energies, XMM- 
COSMOS will provide a valuable sample of absorbed sources 
to test X-ray background (XRB) synthesis model predictions. 
Moreover, to understand the nature of the XRB sources, it is 
very important to have a detailed, cosmic variance free, mea- 
surement of the amplitude of the logN-logS relations in several 
energy bands. It is also worth noting that XMM-COSMOS sam- 
ples with good accuracy the flux range where most of the 
XRB flux is produced (i.e. around S(2-10 keV)~10~ 14 erg 
cm" 2 s" 1 ). Therefore, am ong the medium-deep X-ray surveys 
dBrandt & Hasingerll2005h . XMM-COSMOS has the best com- 
bination of these characteristics to achieve the goals mentioned 
above. 

The XMM-COSMOS survey, with its large area and counting 
statistics, provides a large sample of bright sources where the 
hardness ratio can be measured with good precision. Thanks also 
to the large amount of spectroscopic data in the field it is pos- 
sible to compare, in a reliable way, the optical properties with 
the X-ray properties derived from the hardness ratio analysis for 
large samples of sources. This is particularly important for AGN 
classification into absorbed (Type II) and unabsorbed (Type I). In 
recent years it was realized (Szokoly et al. 2004) that the classifi- 
cations based on optical spectroscopy may be affected by strong 
biases and AGN can be missed or not recognized as such. 
The paper is organized as follows: in Section [2] we present the 
observations and we summarize the data reduction techniques; 
in Section 13.11 we report on the source detection; in Section 
13.21 we present the pointlike source catalog; in Section 13.31 we 
quantify, using Monte Carlo simulations, the completeness of 
the catalogue; in Section [4] we present the logN-logS relations; 
in Section [5] we give an overview of the source content of the 
field using X-ray colour-colour diagrams and the overall results 
are summarized in Section[6] Unless otherwise stated, errors are 
given at the lcr level and we assume a A dominated Universe 
with H =70 km/s/Mpc, Q,„=0.3 and Q A =0.7. 

2. Observations and data reduction 

The XMM-COSMOS survey covers 2.13 deg 2 in the equato- 
rial sky in a region bounded by 9 ;i 57.5'" < a <10' 1 03.5" ! and 
1°27.5' < 6 <2°57.5'. X-ray observations were performed 
during XMM-Newton A03-A04 from December 2003 to June 
2006. The survey consists of a matrix of 5x5 pointings shifted 
by 15' with respect to each other. The matrix of pointings was 
observed in A03 and repeated with a rigid shift of 1' in A04. 
The shift was applied to smooth sensitivity drops introduced by 
the CCD gaps. In Table Q] we present the log of the 55 XMM- 
Newton observations of the COSMOS field. 
Because of charged particle flares, two pointings were com- 
pletely lost, namely 16A and 25 A. The lost times were com- 
pensated for by tuning the exposures in A04. Additionally, two 
pointings (i.e. field 20C and 23C) were re-observed in XMM- 
Newton A06 (May 2007) for 32 ks each to compensate for time 
losses. At the time of writing no further observing campaigns of 
the COSMOS field are planned with XMM-Newton. 



In Paper II we analyzed a first sample of 23 fields observed with 
XMM-Newton during A03 labeled in Table [T] The total expo- 
sure was ~504 ks after the cleaning of the particle background 
flares. The faintest sources in the field have a flux of 7x 10~ 16 erg 
cm' 2 s , 4xl0~ 15 erg crrT 2 s and 9xl0~ 15 erg crrT 2 s in the 
0.5-2 keV, 2-10 keV and 5-10 keV energy bands, respectively, 
while a flux limit of ~1.7xl0~ 15 erg cirT 2 s -1 , ~9.3xl0~ 15 erg 
cm -2 s _1 and ~1.3xl0~ 14 erg crrT 2 s was achieved over 90% 
of the area (1.92 deg 2 ) in the 0.5-2 keV, 2-10 keV and 5-10 keV 
energy band, respectively. The preliminary catalogue based on 
those data consisted of 1390 independent sources and 1281, 784 
and 186 source in the three bands, respectively. We used that cat- 
alog to produce the first XMM-COSMOS logN-logS relations as 
well as the first study of the cosmic or sample variance in X-ray 
surveys. Paper II also contains a detailed section on data analy- 
sis techniques, including event cleaning, image processing, as- 
trometry, source detection and Monte Carlo simulations. In this 
section we briefly summarize the analysis method; we refer the 
reader to Paper II for a detailed description. 
XMM-Newton was operated in imaging mode using the EPIC 
CCD cameras in full frame mode. X-ray event files were 
searched for particle background flares and screened with the 
technique described in Paper EL In order to reduce the instru- 
mental background, the energy channels between 1 .45 ke V and 
1 .54 keV were discarded in both the MOS and PN data. To re- 
move the strong Cu fluorescence features in the PN background 
we also discarded the energy bands 7.2 keV-7.6 keV and 7.8 
keV-8.2 keV. The total scheduled EPIC exposure time was 1464 
ks while, after the background cleaning the sum of the PN good 
time intervals (GTI) was -988 ks and 1207 ks for both MOS1 
and MOS2. 

Due to the slow decrease of the solar activity from its maximum 
in 2000 to its minimum in 2007 dHathawav et aljl9 99), observa- 
tions performed in A03 and in the first part of A04 have a signif- 
icantly higher background level than in the second part of A04 
and the two observations in A06. Event files were processed us- 
ing the XMM-Newton Standard Analysis Software (SAS) ver- 
sion 6.7.0. After the removal of high background intervals we 
searched for and removed hot/dead columns and pixels. Images 
were created in the 0.5-2 keV, 2-8 keV and 4.5-10 keV energy 
bands. In the same bands we created spectral weighted exposure 
maps assuming a power-law model with photon spectral index 
T=2 in the 0.5-2 keV band and T=1.7 in the 2-8 keV and 4.5- 
10 keV bands. 

The 0.5-2 keV exposure map of the XMM-COSMOS survey is 
shown in Fig. Q] while in Fig. [2] we show a false colour X-ray 
image of the entire field. 

In order to compute background maps, we performed a prelim- 
inary source detection using a sliding cell technique. Using a 
threshold of 2.5cr with the XMMSAS software 'eboxdetect', we 
excised all the detected sources from all the images. The result- 
ing images were fitted with a double component model (a flat 
and a vignetted component) to mimic the particle and the X-ray 
sky background. 

Astrometry corrections were estimated as in Paper II by 
cross-correlating highly significant (i.e. det_ml>15, see Sect. 
13.11 ) X-ray sources detected in each pointing, with the cata- 
log of galaxies detected in the I-band by CFHT-MEGACAM 
dMcC racken et al. 2007) and computing the most likely shift us- 
ing the XMM-SAS software "eposcorr". The mean astrometric 
shift is similar to that reported in Paper II, being A(a) ~1.4" and 
A(6) -0.2". 
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Fig. 1. Colour coded vignetting corrected 0.5-2 keV exposure map of the XMM-COSMOS survey. The maximum effective depth 
achieved on the field is ~84 ks (white) and the mean exposure is ~68 ks (green). 



3. Source detection and source catalogue 

3. 1 . Source detection 

We ran a two steps source detection in three energy bands, i.e 
the 0.5-2 keV, 2-8 keV and 4.5-10 keV. 

By using the XMM-SAS tool "eboxdetect" we first ran a sliding 
cell detection to select source candidates in the field. Differently 
from Paper II, we used the 2-8 keV band in place of the 2-4.5 
keV band. This because a reanalysis of the A03 data showed 
that the 2-8 keV band yields a better estimate of the 2-10 keV 
source counts (i.e. in the 2^4.5 keV band we detected 10% fewer 
sources than in the 2-8 keV band). Moreover we determined that 
excluding the 8-10 keV events from the analysis slightly en- 
hanced the signal-to-noise ratio of most the 2-10 keV sources. 
The 4.5-10 keV band has been fully exploited in order to find 
the most absorbed sources. 

If P is the probability that a Poissonian fluctuation of the back- 
ground is detected as a spurious source, the likelihood of the 
detection is then defined as det_ml=- ln(f). All the source 
candidates with det_ml< 4 were discarded. Making use of 
the XMM-SAS tool "emldetect" we then performed a maxi- 
mum likelihood fit of each source candidate to a PSF model 
available in the XMM-Newton libraries. All the sources were 
also fitted with a convolution of a j8-mode l cluster bright- 
ness profile (Cavali ere & Fusc o-Femiariol l 19761) with the XMM- 



Newton PSF, to determine a possible extension in the detected 
signal. A source is classified as extended if the likelihood for the 
y6-model fit exceeds that of the pointlike case of 10 in det_mZ. 
The sources are fitted simultaneously in the 0.5-2 keV, 2-8 keV 
and 4.5-10 keV energy bands and the free parameters of the fit 
are position, flux and extension. 

Moreover, the calculation of the positional uncertainties in each 
band also makes use of the information in other bands and thus 
source positioning is extremely accurate in all the energy bands 
( see also discussion in Paper II). In Fig. [3] we show the dis- 
tribution of the statistical uncertainties on the source positions 
in arcsec as output by the emldetect task. The median statis- 
tical astrometri c uncertainty, inc luding also a systematic error 
of 0.75" (see iBrusa et al.ll2009l for a detailed discussion), is 
1 .77"B The reliability of the estimated source positions is con- 
firmed by the distribution of the offset between X-ray sources 
and optical counterparts. Count rates estimated in the 2-8 keV 
and 4.5-10 keV energy bands were extrapolated into 2-10 keV 
and 5-10 keV fluxes, respectively. In these bands we com- 
puted energy conversion factors (ECF) by assuming a power- 
law spectrum with spectral index T=1.7 and Galactic column 
density N H =2.5xl0 20 cm 4 . In the 0.5-2 keV band, we directly 



1 Similar results also have been obtainedvia Monte Carlo simula- 
tions. 
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Fig. 2. False colour X-ray image of the COSMOS field: red,green and blue colours represent the 0.5-2 keV, 2^1.5 keV and 4.5-10 
keV energy bands, respectively. 



converted the count-rate into fluxes assuming a spectral index 
T=2.0 and Galactic column density N ff =2.5xl0 20 cm" 2 . The 
choice of these spectral in dices is driven by the findings of 
Main ieri et al.l d2007l 12008). They measured an average spec- 
tral index < T >~1.7 for the XMM-COSMOS sources. In the 
soft band we have chosen a steeper index to take into account 
the contribution of the soft excess. More over these values of T 
are widely used in the literature (see e.g. lHasinger et al1 [l993; 
iBaldi et al]|2002l) and therefore this choice has also the scope of 
a better comparison with previous works, especially when com- 
paring the logN-logS relations. 

The adopted ECFf§ are 10.45 cts s" 1 / 10" 11 erg cm" 2 s" 1 , 2.06 
cts s" 1 / 10" 11 erg cm" 2 s" 1 and 1.21 cts s" 1 / 10" 11 erg cm" 2 
s" 1 in the 0.5-2 keV, 2-8 keV and 4.5-10 keV energy bands, 
respectively. All the sources with a maximum likelihood param- 
eter det_ml>10 in at least one band have been included in the 
present catalog. This threshold corresponds to a fraction of ex- 
pected spurious sources of the order of 1.5% in the 0.5-2 keV 
band and ~0.5% in the other energy bands. Since in this work 
we used a more conservative detection threshold than in Paper 
II (det_ml>6), the fraction of spurious sources has been signifi- 
cantly reduced. 



Significant detections have been achieved only in a subset of the 
energy bands. In the bands where the detection is not signifi- 
cant, we comp uted 1 crupp er limits of the counts using the pre- 
scriptions of Narskv (2000). Given M counts actually measured 
in a region of 30" at the position of the source and B back- 
ground counts (estimated from our background maps), the la- 
upper limit is defined as the number of counts X that gives the 
probability of observing M (or fewer) counts equal to the formal 
68.3% Gaussian probability: 

P(< M,X + B) = P Cfl „, 5 (68.3%). (1) 
Assuming Poissonian statistics, this equation becomes: 

Poaus^e-^^^^.. (2) 

!=0 l ' 

By solving Eq.|2]iteratively in the case of Pcau.ss-0.683, we ob- 
tained the lcr upper limit X. The upper limits were then con- 
verted into count-rates and fluxes by diving by the exposure map 
and then applying the ECFs. We removed from the catalogue 
about 20 sources lying close to clear artifacts of the image (i.e. 
field and pointing boundaries or unremoved hot pixels). With the 



2 The ECF values also take into account the energy channels dis- 
carded to decrease the background. 



3 We checked that a 30" aperture gives the best agreement between 
aperture photometry and the maximum likelihood PSF fitting technique. 
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1 42 36.0 


31.9 


7.0 


9.2 


20B 


1178 


0302352001 


2006-05-15T17:24:58 


09 59 22.4 


1 42 36.0 


19.9 


4.9 


16.6 


20C* 


1356 


0501170101 


2007-05-06T00:23:52 


09 59 22.4 


1 42 36.0 


33.9 


32.0 


33.3 


21A" 


916 


0203362101 


2004-12-09T07:16:01 


09 58 26.4 


2 42 36.0 


62.6 


60.3 


61.7 


22A« 


898 


0203362201 


2004-1 1-03T06:02:44 


09 58 26.4 


2 27 36.0 


30.9 


28.0 


30.5 


22B 


1176 


0302352201 


2006-05-1 1T19:47:08 


09 58 30.4 


2 27 36.0 


21.9 


7.0 


10.8 


23A 11 


992 


0203362301 


2005-05-09T19:01:30 


09 58 26.4 


2 12 36.0 


30.9 


1.3 


28.1 


23B 


1176 


0302352301 


2006-05-12T02:30:29 


09 58 26.4 


2 13 36.0 


21.9 


4.3 


7.4 


23C* 


1362 


0501170201 


2007-05- 18T03: 17:39 


09 58 26.4 


2 13 36.0 


36.0 


28.1 


33.9 


24A« 


992 


0203362401 


2005-05-10T04:14:50 


09 58 26.4 


1 57 36.0 


30.9 


17.4 


23.0 


24B 


1175 


0302352401 


2006-05-09T19:36:30 


09 58 30.4 


1 57 36.0 


24.9 


0.2 


21.9 


24C 


1190 


0302353201 


2006-06-09T01:36:12 


09 58 30.4 


1 57 36.0 


19.3 


9.7 


14.6 


25A" 


992 


0203362501 


2005-05-10T13:28:ll 


09 58 26.4 


2 42 36.0 


31.9 


0.0 


0.0 


25B 


1175 


0302352501 


2006-05-10T03:09:54 


09 58 26.4 


1 43 36.0 


24.5 


22.6 


23.9 


25C 


1190 


0302353301 


2006-06-09T07:36:ll 


09 58 26.4 


1 43 36.0 


18.9 


11.6 


14.3 



Table 1. The XMM-Newton observation log of the XMM-COSMOS survey. From left to right: Field ID,revolution, OBSJD, 
observation start, right ascension, declination, duration of the exposure, Good Time Interval (GTI) for the PN and MOS camera, 
respectively. 

Fields observed in XMM-Newton A03 presented by Hasinger et al. 1 2007) and used for Paper II. 
b : Fields observed in XMM-Newton A06. 



Band 


Total detections 


Single-Band detections 


S/im 

erg cm~ 2 s~'/10~ 15 


S50% 

erg cnT 2 s~'/10 15 


Sso% 

erg cnT 2 s~'/10~ 15 


erg cnT 2 s~'/10~' 5 


0.5-2 keV 


1621 


771 


0.50 


1.00 


1.70 


3.00 


2-10 keV 


1111 


237 


2.50 


5.60 


9.30 


15.00 


5-10 keV 


251 


5 


5.10 


11.00 


13.0 


20.00 



Table 2. Summary of the total detections, single band detections, faintest flux limits, flux limits at 50%, 90% of the total area and 
flux limits observable on the full area in the 0.5-2 keV, 2-10 keV and 5-10 keV energy bands, respectively. 
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12 3 4 

Astrometric error (arcsec) 

Fig. 3. The distribution of the positional uncertainties of the 
XMM-COSMOS detections. 



method described above we selected a total of 1887 independent 
sources. Each source has been named with a unique ID number. 
1621 sources have been detected in the 0.5-2 keV energy band, 
while 1111 and 251 sources are detected in the 2-10 keV and 
5-10 keV band, respectively. The number of sources with a sig- 
nificant detection in only one band is 771 for 0.5-2 keV band, 
237 and 5 for the 2-10 keV and 5-10 keV bands, respectively. 
The faintest sources in the field have fluxes of 5.0xl0~ 16 erg 
cm -2 s , 2.5xl0~ 15 erg crrT 2 s~'and 5.1xl0~ 15 erg cm -2 s _1 in 
the three energy bands. 

A summary of the source detection results is shown in Table [2] 
Thanks to our PSF fitting technique we were able to detect 
109 additional extended sources. The catalog of the extended 
sources, together with a detailed and more extensive analysis 
of their properties will b e presented in a forthcoming paper by 
Finog uenov et all d2008l) . 



3.2. Source catalogue 

In Table. |3] we show, as an example, the first 50 entries of the 
catalogue as they appear on-line. The table is structured as 
follows: 

Column 1: IAU Name, Column 2: XID, Column 3: a (deg), 
Column 4: 6 (deg), Column 5: Positional error (arcsec), 
Column 6: 0.5-2 keV flux (erg cm" 2 s^/lO" 14 ), Column 7: 0.5- 
2 keV net counts , Column 8: 0.5-2 keV likelihood parameter 
det_ml, Column 9: 0.5-2 keV background counts (cts/pixfl 
Column 10: 0.5-2 keV vignetting corrected exposure (ks), 
Column 11: 2-10 keV flux (erg cnr 2 s^/lCT 14 ), Column 12: 
2-8 keV net counts , Column 13: 2-8 keV likelihood parameter 
det_ml, Column 14: 2-8 keV background counts (cts/pix), 
Column 15: 2-8 keV vignetting corrected exposure (ks), 
Column 16: 5-10 keV flux, (erg cirT 2 s^'/lCT 14 ), Column 17: 
4.5-10 keV net counts, Column 18: 4.5-10 keV likelihood 



parameter det_ml, Column 19: 4.5-10 keV background counts, 
Column 20: 5-10 keV vignetting corrected exposure (ks). 
The interactive and machine readable full 
version of the catalog can downloaded at 
http : //irsa . ipac . caltech . edu/data/COSMOS/. For 
sources with no significant detection in a band, we list upper 
limits with negative values of flux. In this case, we also quote a 
value of cts=0, det_ml=-l and Bkg=-1 (background) in the band 
where the detection is not significant. 

The flux errors are the statistical uncertainties estimated from 
the maximum likelihood and do not include uncertainties 
introduced by the choice of the spectral model to estimate the 
flux. We determined that by varying by AT=0.3 the spectral 
index assumed in computing the fluxes, the resulting variation 
of the flux estimate is of the order of 2%, 9%, and 4% in the 
three bands. 

The Chandra coverage of the inner area of the XMM-COSMOS 
field (El vis et al.l l2009) offers a unique possibility to investigate 
the effect of source confusion in our catalog. 
The Chandra field covers about half of the XMM-COSMOS 
field. Of the 1887 XMM-sources with det_ml>10, 946 (50.1%) 
have been observed by Chandra with an exposure longer than 
30 ks, and 876 of them are prese nt in the C-COSMOS point-like 
source catalog dElvis et al.ll2009h . Twenty-four of the 876 XMM 
pointlike sources with Chandra coverage (2.7%) are actually 
resolved into two different Chandra sources, which lie between 
2 and 10 arcsec from each other and have been blurred by the 
XMM large PSF. We then used the Chandra source counterpart 
positions of these 24 "blended sources" (48 different positions) 
as the input catalog for emldetect and we fitted these sources 
keeping the position parameter fixed at the Chandra value. As a 
result only 2/24 XMM-COSMOS sources have been deblended 
into 4 XMM-Newton sources, namely XID #67 and XID #82, 
while the remaining 22 sources have been detected again as a 
single XMM-COSMOS source with properties consistent with 
those presented in the catalogufl Therefore we can conclude 
that our sample contains <2.7% of the sources which could 
be resolved into two sources at the Chandra-COSMOS flux 
limit. In Fig. [4] we show the XMM-Newton image of a region 



The pixel scale is 4"/pix 



5 We kept these sources as single entries in the catalog for self con- 
sistency with our statistical analysis 
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Fig. 5. Left Panel : The sky coverage versus flux for the 0.5-2 keV, 2-10 keV and 5-10 keV bands, represented by red, green and 
blue solid lines. The horizontal solid and dashed lines show the 90% and 50% completeness levels. Right Panel : The 0.5-2 keV 
sensitivity map of XMM-COSMOS in erg cirT 2 s . The map is plotted in colour coded scale from 5xl0~ 16 erg cirT 2 s _I (magenta) 
to 3xl0~ 15 erg cirT 2 s _1 (red). 



containing the two deblended sources, which by chance are 
close to each other. 

3.3. Monte Carlo simulations, sky coverage and sensitivity 
maps 

In order to estimate the sky coverage of our survey, we per- 
formed Monte Carlo simulations as described in Paper II. The 
precision of the photometry as well as positional accuracy were 
also discussed in Paper II. Here we give an overview of the pro- 
cedure adopted for the production of random X-ray sky images 
and their analysis. 

Twenty series of 55 XMM-Newton images were created with 
the same pattern, exposure maps and background levels as the 
real data. We produced 20 random input source catalogs with 
sources randomly placed in the field of view and fluxes dis- 
tribute d according to the AGN logN-logS distributions predicted 
by the iGilli et al.l d2007l) XRB synthesis model. The input fluxes 
were converted into count-rates by folding through the response 
matrices the same spectral model assumed to compute fluxes and 
to weight exposure maps. The counts of the sources were then 
convolved with XMM-Newton PSF templates available in the 
XMM-Newton calibration database and reproduced on the de- 
tector. We then applied to the simulated fields the same source 
detection procedure used in the real data producing 20 indepen- 
dent output catalogs. 

The sky coverage is then obtained by dividing the number of 
detected sources at each flux by the number of input sources 
an d rescaling for a to tal area of 2.13 deg 2 . By using as a model 
the lGilli et al](l2007l) logN-logS, it is possible that the simulated 
logN-logS could be slightly different from the real logN-logS. 
This could introduce some biases in the es t imatio n of the ef- 
fective area. However, Schm itt~& Maccacarold 19861) showed that 



the effect of a different slope of the logN-logS is negligible when 
the threshold of the source detection is higher than 3-4cr. 
The sky coverage in the three energy bands under investigation 
is plotted in the left panel of Fig.|5]As a result of the simulations, 
we obtained that 90% of the survey area is sensitive to flux lim- 
its of ~1.7xl0~ 15 erg cirT 2 s -1 , ~9.3xl0~ 15 erg cirT 2 s -1 and 
~1.3xl0~ 14 erg crrT 2 s _1 in the three energy bands, respectively. 
Additionally, we determined that the survey is sensitive, over the 
full field of view (i.e. 2.13 deg 2 ), to fluxes of ~3.0xl0~ 15 erg 
cm -2 s _1 , ~1.5xl0~ 14 erg cirT 2 s _1 and ~2.0xl0~ 14 erg cm -2 
s _1 , respectively. As mentioned above, the fluxes of the input 
spectrum are converted into count-rates by assuming a single 
spectrum for all the sources. This could in principle bias the es- 
timates of the sensitivity limits. In order to test the effect of a 
variation of the mean spectral index in the estimate of the sky 
coverage, we changed the spectral indices by AT = +0.3. In this 
way the estimate of the flux limit changed by <2% in the soft 
band, while this variation was of the order of 9% and 4% in the 
2-10 keV and 5-10 keV band. The sky coverage is thus almost 
insensitive to a change of the spectral shape in the 0.5-2 keV 
and in the 5-10 keV band. In the 2-10 keV we estimated that a 
9% uncertainty in the flux limit could introduce an overall uncer- 
tainty of ~5% in the logN-logS. Such an uncertainty is however 
smaller than the typical uncertainty on the source counts. 
In order to map the sensitivity across the field of view we pro- 
duced sensitivity maps of the XMM-COSMOS survey in all the 
energy bands by reversing our source detection analysis. By us- 
ing our estimated background maps and exposure maps we eval- 
uated, according to the Poisson statistic, the minimum number 
of counts necessary to have a detection with detjnl >10. The 
number of counts have been evaluated in cells of 3x3 pixels and 
corrected for the fraction of the PSF falling out of the cell. The 
resulting count-limit maps have been divided by the exposure 
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maps and converted into flux limit maps using the ECF. 
As an example, the resulting 0.5-2 keV band sensitivity map 
is plotted in colour scale in the right panel of Fig. The map 
is in excellent agreement with the sky coverage plot obtained 
via Monte Carlo simulations. As one can notice, almost all the 
central area (~1.8 deg 2 ) has a quite homogeneous flux limit 
~1.7xl0~ 15 erg cirT 2 s _1 . The northern central part of the field 
shows an area of ~0.5 deg 2 , having a flux limit of the order of 
~8xl0~ 16 erg cirT 2 s . It is worth noting that the deepest part 
of the field is located in the northeastern part of the field and the 
flux limit of ~5.0xl0~ 16 erg cm -2 s _1 in the sensitivity maps is in 
agreement with the predictions of the Monte Carlo simulations 
and the output of the source detection. 



4. logN-logS relations 

Using the sky coverage we produced the cumulative logN-logS 
relations in the three energy bands under investigation by using: 



m>s) 



N s . 

5^ 



,-2 



(3) 



where N(> S ) is the total number of detected sources in the field 
with fluxes greater than S and £2, is the sky coverage associated 
with the flux of the i' A source. The variance of the source number 
counts is therefore defined as: 



N s , . 



(4) 



The cumulative number counts, normalized to the Euclidean 
slope (i.e. multiplied by S 15 ), are shown in Figure [6] [7] and [8] 
in the 0.5-2 keV, 2-10 keV and 5-10 keV energy ranges, re- 
spectively. The logN-logS relations are also presented in Table 
|4] From left to right: Flux, Number-counts and area in the 0.5-2 
keV, 2-10 keV and 5-10 keV energy band, respectively. 

In order to parametrize our relations, we performed a max- 
imum likelihood fit to the unbinned differential counts. We as- 
sumed a broken power-law model for the 0.5-2 keV and 2-10 
keV bands: 



dN 

n(S) = — = 
as 



AS' 01 S >S b 
BS" 2 S < S b , 



where A = BS? 1 '" 2 is the normalization, a\ is the bright end 
slope, ff2 the faint end slope, Sb the break flux, and S the flux 
in units of 10~ 14 erg cirT 2 s -1 . Notice that, using the maximum 
likelihood method, the fit is not dependent on the data binning 
and therefore we are using the whole dataset. Moreover, the nor- 
malization A is not a parameter of the fit, but is obtained by im- 
posing the condition that the number of expected sources from 
the best fit model is equal to the total observed number. 

In the 0.5-2 keV energy band the best fit parameters 
are an =2.40+0.05, ar 2 =1.60!^, S^l.OO^'xlO^erg cnr 2 
s _1 and A=141. These values are consistent with those mea- 
sured in paper II while the normalization is lower than the value 
(A=198) derived in paper 10 However, with this fitting method 
the normalization is not a fit parameter and it is strongly depen- 
dent on the best fit values of the bright end slope and on the cut- 
off flux. One can indeed notice that the best fit values of the a\ 
and Si, parameters are somewhat changed with respect to paper 



6 Note that in paper II we gave the normalization of the cumulative 
distributions 
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Fig. 6. The 0.5-2 keV logN-logS of the XMM-COSMOS 
(red dots) so urces compared with the ROSAT medium sensi- 
tivity survey {Hasinger et al. 1993, blue dot das hed line), com- 
bined ROSAT, XMM-Newton, Chandra sources (lHasinger et al 



2005, green dashed line), the 2Ms CDFS ( 1 <x error tie. lLuo et al 
2008, magenta contin uous line), the 2 Ms CDFN (l<x error 
tie, iBauer et all 12004 p ink dot - dashed line), the XMM- 
Newto n Lockman hole ( Brunner et"aT]|2008l blue circles), the 



AXIS (ICarrera et al.l2007l cyan triangles), the HELLAS2XMM 
(Bald i et al.l2002l black pen tagons) and the extended CDFS (l<x 
error tie. lLehmer et al.l200 5, black continuous line) surveys. The 
source number counts are plotted multiplied by (S/10 14 ) 15 in or- 
der to highlight the deviations from the Euclidean behavior. 



II. The bright end slope varied from 2.6 in paper II to 2.4 in the 
present work and the cut-off flux varied from ~1.55xl0" 14 erg 
(*) cirT 2 s _1 to 1.00xl0~ 14 erg cirT 2 s _1 . However, a comparison of 
the amplitude of the source surface density measured in paper 
II with that measured here can be performed if we measure the 
model predicted source counts at fluxes fainter than the knee. If 
we take 2xl0~ 15 erg cirT 2 s _1 as a reference flux, in paper II we 
had a source surface density of 478 deg~ 2 while here we measure 
479 deg -2 . We can therefore conclude that the 0.5-2 keV logN- 
logS obtained in paper II and in this work are in good agreement. 
In the 2-10 keV band the best fit parameters are a\ =2.46+0.08, 
0-2=1-55+0.18, 5i,=1.05±0.16 10" 14 erg cnr 2 s 1 and A=413. 
Since the best fit parameters are similar to those of Paper II, we 
can directly compare the normalizations of the logN-logS. The 
normalization derived in this work is 10% higher than that mea- 
sured in paper II. This effect is partly due to the sources missing 
in the 2^1.5 keV band and detected in the 2-8 keV which were 
not considered in the analysis of paper II. Moreover, extrapola- 
tion of the 2-4.5 keV count-rate into the broader 2-10 keV band 
is more affected by uncertainties on the true source spectral slope 
and provides wrong count-rate estimates especially for the most 
absorbed sources. 

In the 5-10 keV energy band we did not find any signif- 
icant break in the slope. We therefore fitted the data using 
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Fig. 7. The 2-10 keV logN-logS of the XMM-COSMOS 
(red dots) sources compared with the combined Chandra , 
XMM-Newton and ASCA sources (iMoretti et all I2003L 
blue dashed line), the HELLAS BeppoSAX dGiommi et alj 
2000, black h exagons) the 2 Ms CDFS (lcr error tie, 



Luo et al]|2008l m asenta continuous line), the HE LLAS2XMM 
(Bald i et all 120021 green pentagons), the AXIS (ICarrera et alj 
2007, cyan triang les), the extended CDFS ( lcr error tie, 
Lehmer et al. 2005, black continuous line) and the Lockman 
hole dBrunner et alj 120081 blue open circles) surveys. The 
black - dashed - line are the lcr confi dence contours of th e best 
fit to the logN-logS of the ASCA data dCagnoni et al.ll998l) . The 
source number counts are plotted multiplied by (S/10 14 ) 15 in or- 
der to highlight the deviations from the Euclidean behavior. 



a single power-law in the form of n(S)=AS a> and obtained 
ai =2.38+0.05 and A=130. 

In the 5-10 keV both the normalization and the slope are consis- 
tent within lcr with the values obtained in paper II. 



4. 1 . Comparison with previous surveys 

In Fig. [6] I7I8I we compare our logN-logS with the results of 
previous surveys. A visual inspection of the data shows that 
the XMM-COSMOS source counts are in general agreement, 
within lcr, with all the previous measurements. In the 0.5-2 
keV band source counts of all the surveys agree with our 
measurements, with the only exception of the bright end of the 
Lockman Hole logN-logS. The reason for such a discrepancy 
is that the location of the Lockman Hole survey was chosen 
on purpose near a concentration of bright sources to improve 
the accuracy of the ROSAT star tracker in order to achieve a 
better astrometry (G. Hasinger, private communication). This 
had the result of artificially increasing the source counts at the 
bright end of the relation. The comparison with other surveys 
is consistent with the error bars and with the counts in cell 
fluctuations predicted in Paper II. We a lso co mpared our results 
with the recent work of Mat eos et alj d2008l) who performed a 



Fig. 8. The 5-10 keV logN-logS of the XMM-COSMOS 
(red dots) sourc es compared with the HELLAS2XMM 
dBaldi et all 120021 green pe ntagons), the 1 Ms CDFS (lcr 
error tie, iRosati et ail 120021 masenta continuous line), the 
HELLAS-B eppoSAX dFiore et al.N200lL black hexagons), the 
ELAIS SI dPuccetti et alj 120061 blue stars), the XMM-HBS 
(Delia Ce ca et al] 12004, blue triangles) and the Lockman Hole 
dBrunner et alj 120081 blue open circles), surveys. The source 
number counts are plotted multiplied by (S/10 14 ) 1 5 in order to 
highlight deviations from the Euclidean behavior. 



detailed analysis of the logN-logS of X-ray sources detected in 
1129 XMM-Newton archival observations. By com paring the 
data o f Table [4] with those shown in Table 3 of Mateos et al. 
(2008) we found lcr agreement in almost all the data bins. 
In paper II we showed that the fluctuations of the source counts 
are proportional to the actual number of sources in the field and 
to the amplitude of the angular auto-correlation function of the 
X-ray sources. Therefore, assuming a universal shape of the 
autocorrelation function, we expect that the surveys showing 
the largest deviations from the mean value of the source density 
are the pencil beam surveys (i.e. area <0.2 deg 2 ) at their bright 
end. Moreover, with XMM-COSMOS , fluctuations introduced 
in previous shallow surveys by low counting statistics and 
by random sampling of a few large structures in pencil beam 
surveys are largely suppressed. With the same formulas used 
in Paper II, we estimate the effect of the cosmic variance to be 
<5% on the normalization of the XMM-COSMOS logN-logS 
and that the new data do not change the results shown in Paper 
II. 

Also in the 2-10 keV energy bands we do not note any signifi- 
cant deviation from previous works with the exception only of 
the Lockman Hole and the two faintest bins of the AXIS counts. 
A lso in this band our d ata are in good agreement with the results 
of lMateos et ai1d2008l) . 

Fig. [8] suggests that the fluctuations of the source counts in 
the 5-10 keV band are much larger than in the other bands. 
This is due to the fact that as discussed above and in paper II, 
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when we deal with low source surface density, the impact of 
the sample variance becomes significantly high. However in 
this band our data are statistically consistent with most of the 
data from other surveys. Also in this energy band, the deviation 
of the COSMOS data from those of the Lockman Hole is due 
the higher number of bright sources in that particular field. Our 
source counts are 1 0-1 5 % h igher than those of the ELAIS-S1 
survey dPuccetti et al.l 12006). As an example, using Eq. 10 of 
paper II, we determine that at the faint end of the 5-10 keV 
band, fluctuations due to the cosmic variance are of the order 
of the 20-40 percent, depending on the survey size. At the 
bright end large deviations of more than a factor of two are 
still allowed by the sample variance. This is also visible in 
Fig. [8] where at the bright end the Beppo-SAX counts exceed 
the XMM-COSMOS counts by about a factor of two though 
remaining statisti cally consistent with each other. 
iKim et all (120071) reported the results of a broken power law 
fit to the logN-logS from different surveys available in the 
literature. They also reported measurements of the CHAMP 
survey which is a compilation of Chandra archival data for a 
total sky coverage of 9.6 deg 2 , with a depth about one order of 
magnitude fainter than XMM-COSMOS. On average the bright 
end slopes are consistent with a Euclidean rise in all the surveys. 
The faint end slopes are of the order of ~ 1.5- 1.6 in the 0.5-2 
keV band and span from ~1.3 to a<x ~2.0 with a mean of 
02 ~1. 6-1.7. A larger spread is reported for the cut off fluxes. 
Although the spread in this parameter is quite large, our data are 
consistent with the average values reported in the literature for 
this parameter. 



4.2. Extragalactic X-ray source number counts and 
comparison with models 



0.3 





io- 12 




Fig. 10. Upper panel The ratio between the iGilli et al.l (l2007h 
model logN-logS relations to the observed sour ce counts in 
XMM-COSMOS and in the Chandra deep fields dRosati et all 
2002; iBauer et all 12004 . From bottom to top in the 0.5- 
2 keV, 2-10 keV and 5-10 keV energy bands. The XMM- 
COSMOS datapoints are plotted in red, while in black and 
blue we plot the CDFN and CDFS , respectively. Bottom p anel: 
Same as U pper panelbut using the lTreister & Urry (2006) XRB 
model. 



Fig. 9. The 0.5-2 keV flux distribution of sources classified as 
AGN or extragalactic (black) and stars (red). 



We used our logN-logS relations to test the most recent ex- 
tragalactic XRB synthesis models. In order to compare our data 
with the XRB m odel, we est i mate t he fraction of sources clas- 
sified as stars by iBrusa et al.l (120091) . In the 0.5-2 keV band we 
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identified 74/1621 (i.e. ~4.5%) sources classified as stars, while 
these are 17/1111 (i.e. -1.5%) and 3/251 (i.e. -1.1%) in the 2- 
10 keV and 5-10 keV bands, respectively. In figure[9]we plot the 
normalized distributions of the fluxes of stars and extragalactic 
sources in the 0.5-2 keV band. Since the two distributions are 
similar we can conclude that stars in the XMM-COSMOS flux 
range affect the 0.5-2 keV logN-l ogS only by increas ing the ex- 
tragalactic source counts by ~5%. lMateos et al] (12008) measured 
a flux dependent fraction of stars, with higher fractions than ours 
at bright fluxes where XMM-COSMOS is undersampled. By ex- 
cluding the source classified as stars, we derived the logN-logS 
relations for extragalactic sources only. 

In the upper panel of Fig. [10] we plot the ratio of the XMM- 
COSMOS logN-logS relations t o the predi ct ions o f the XRB 
population synthesis model of iGilli et ail d2007l) (hereafter 
model I) , while in the bottom pan el we plot for compari son the 
ratio of the data to the model of iTreister &Urrvl d2006l) (here- 
after model II). In both models the XRB spectrum is dominated 
by obscured AGN which outnumber unobscured ones by a factor 
3-4 at low X-ray luminosities (logL^ < 44). The cosmological 
evolution is similar and parametrized using the most recent de- 
terminations of the AGN luminosity function (Ueda et al. 2003; 
La Franca et al. 2005; Hasinger et al. 2005). In both models the 
obscured fraction decreases towards high luminosity. The lumi- 
nosity dependence is stronger in Treister et al. (2006) who also 
allow the obscured fraction to increase at high redshifts. The ab- 
sorption distribution is peaked around \ogNn ~ 23.5 in Gilli 
et al. (2007), while it remains rather flat above \ogNn - 22 in 
Treister et al. (2006). They also differ in the adopted XRB in- 
tensity around the 30 keV peak. The Gilli et al. (2007) model 
is tuned to fit the HEAO-1 level, consistent, within 10%, with 
recent BeppoSAX (Frontera et al. 2007) and Swift BAT (Ajello 
et al. 2008) measurements, while Treister et al. renormalize the 
HEAO-1 intensity upward by a factor 1.4 to better match the 
extrapolation of lower energy (< 10 keV) data (i.e. De Luca 
& Molendi 2004). Moreover, for thi s paper we adopt a modi- 
fied version of the IGilli et ail d2007l) modeQ which takes into 
account t he decline o f the sp ace density of AGN at z>3 dis- 
cussed by Brusa et al. (2008). In order to test the models over 
a wider range of fl uxes we als o plotted the da ta of the CDFN 
dBauer et al]|2004 and CDFS dLuo et al.ll2008t) surveys. By re- 
stricting our analysis to fluxes larger than 10~ 16 erg cm" 2 s" 1 , the 
contr ibution of normal galaxy counts is negligible ( Ranall fet al] 
120031) . The results of this comparison can be summarized as fol- 
lows: 

- In the 5-10 keV energy band, both models reproduce well 
the XMM-COSMOS logN-logS, while the CDFS counts 
show a systematically different slope from that of the pre- 
dicted relation. However, because of the small effective area 
of Chandra above 5 keV (i.e. -200 cm 2 at 6.4 keV), the 5-10 
keV CDFS logN-logS may suffer from significant systematic 
uncertainties . 

- In the 2-10 keV energy band, the models repro- 
duce quite accurately the XMM-COSMOS data, although 
model II slightly (i.e. - 10%) overpredicts the XMM- 
COSMOS counts. The CDFN counts show a systematically 
higher normalization than those of the models (up to 40% at 
faint fluxes) and of the COSMOS and CDFS data. 



7 The predictions of the model can be retrieved on line 
at http://www.oabo.inaf.it/~gilli/counts.html using the POMPA 
COUNTS software (POrtable Multi Purpose Application for the AGN 
COUNTS). 



- In the 0.5-2 keV band both models show significant devi- 
ations from both data sets. Source counts estimated from 
model I show a systematically steeper slope than the data. 
On average, model I deviates from the observations by about 
10-15% in the flux range 10~ 16 erg crrT 2 s^'-lO -13 erg cm" 2 
s _1 . Model II, on the other hand, systematically underesti- 
mates the source counts while a visual inspection shows a 
good agreement with the observed slopes. In this case de- 
viations between the data and the model are of the order of 
20% at the XMM-COSMOS fluxes, while at fainter fluxes 
the deviations are larger and of the order of 30%-40%. Both 
models underestimate the observed counts by -30% at fluxes 
greater than -10~ 14 erg cirT 2 s _1 

In summary, the hard X-ray observations are very well repro- 
duced by both models with an accuracy of -10%. In the soft 
band the agreement between the predicted and the observed re- 
lations is not as good as in the harder energy bandsS . The level 
of the discrepancy, however, is small (-20%) and such that it can 
be easily accommodated by slight variations of the XRB model 
parameters. This band is in fact more sensitive to the effect of 
absorption and therefore a fine tuning of absorption in AGN is 
required in the models. Moreover, this ban d contains a larger 
fraction of high-z objects dBrusa et all 2009). Therefore, the fact 
that the two models assume a somewhat different absorption evo- 
lution and XRB spectrum can, in the first instance, explain the 
different source count predictions. We can conclude that at the 
flux limits of the XMM-COSMOS survey, XRB synthesis mod- 
els can reproduce the observations with a precision of 10%-20%. 



5. 



X-ray colours of the X-ray sources 



The X-ray colours or hardness ratios are defined as 

Bj — B] B'i — Bo 

HRi = — and HR 2 



B 2 + Bi 



B 3 +B 2 



(6) 



where B\, B2, and B3 refer to the vignetting-corrected count rates 
in the 0.5-2 keV, 2-10 keV and 5-10 keV bands, respectively. 
By construction, both HRI and HR2 can assume values between 
-1 and 1. 

Fig.QT]displays the HR1-HR2 plot of 212 sources for which the 
lcr error on both HRI and HR2 is <0.25 and for which a high 
quality optical spectrum is available. The plot also contains a 
grid of the expected values of HR 1 and HR2 for different spectral 
models. In particular we considered a simple power law model 
with a spectral index in the interval T=0-h3 and with a column 

de nsity log(N» ) =0-H23 cm' 2 

In iBrusa et all d2009l) . and iTrump et al] (|2009), extragalactic 
sources are classified into 4 main categories: 

- Type I AGN, if the optical spectrum shows evidence of broad 
(FWHM> 2000 Km s" 1 ) emission lines; 

- Type II AGN, if the optical spectrum shows evidence of nar- 
row, high-ionization emission lines and/or AGN diagnostic 
diagrams; 

- Emission line galaxy, if the optical spectrum is dominated by 
a galaxy continuum plus emission lines but without secure 
AGN indicators; 

- Absorption line galaxy if the optical spectrum is dominated 
by a galaxy continuum plus absorption lines; 



8 Note that the inclusion of a decline in the space density of AGN 
at high-z in model I affects mostly the 0.5-2 keV energy bands. In the 
harder bands the predicted number counts are comparable with or with- 
out a high-z space density decline. 
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Fig. 11. X-ray colour-colour diagram in the XMM- 
COSMOS survey. Colours are denned in the text. XID #2608 
has been plotted with its error which also represents the typical 
amplitude of the uncertainties in the plot. The grid represents the 
places in the HR1-HR2 plane of sources with single power-law 
spectra with T=0h-3 with absorption with a column density 
log(N//)=0H-23 crrT 2 . The green line marks the region occupied 
by candidate Compton thick-AGN [i.e. log(N// >24) cirT 2 ] and 
the marks on top of it represent 1%, 3% 10% and 30% level of 
leaking flux, from top to bottom. We represent Type I AGN, 
Type II AGN, emission line and absorption line galaxies as 
blue filled circles, red empty triangles, cyan empty exagons, 
green filled squares, respectively. 



Details of the optic al classificat i on of X-ray sources are exten- 
sively discussed in Brus a et al. (2009), therefore we limit our 
analysis to the X-ray properties of these sources. 140/212 are 
classified as Type I, 32/212 as Type II, 30/212 as emission line 
galaxies and 10/212 as absorption line galaxies. Note that with 
the exception only of 7 objects, all the sources have an estimated 
2-10 keV X-ray luminosity log(Lx)>42 erg/s, with most of them 
having log(Lx)>43 erg/s (see Fig.[T2l. The adopted cuts on the 
errors on the HR preferentially select unabsorbed to moderate 
absorbed AGN, biasing the sample against normal galaxies, star- 
forming galaxies and the most obscured AGN. 
Type I AGN {blue filled circles) cluster in a region around 
HRl=-0.5 and HR2=-0.5 with a relatively small dispersion, cor- 
responding to a typical X-ray spectrum dominated by a power- 
law continuum with very low absorption. Only a few Type 
I sources have X-ray colours typical of Type II sources (i.e. 
HR> -0.1 which corresponds to N# >10 22 cm 2 ). This fraction 
(~2%) is consistent with the results from X-ray spectral analy- 
sis on a subsample of XMM-COSMOS sources (Mainieri et al. 
2007) but at variance with previous works on the fraction of X- 
ray absorbed Type I AGN at comparable X-ray luminosity (see 
e.g. Brusa et al. 2003, Perola et al. 2004, Page et al. 2004) which 
reported values as large as 10%. However such a low fraction 
may be a consequence of the selection effect mentioned above. 



Fig. 12. The 2-10 keV X-ray luminosity vs. HR1 for the sources 
fullfilling the HR error selection. We represent Type I AGN, 
Type II AGN, emission line and absorption line galaxies as 
blue filled circles, red empty triangles, cyan empty exagons, 
green filled squares, respectively. 



On the other hand, type II AGN (red empty triangles) fill most 
of the HR1 -range, corresponding to observed frame absorption 
up to 10 23 cm 2 . The fraction of Type II AGN with X-ray colours 
typical of Type I AGN (HR1< -0.3) is -30%. This is consistent 
with the fraction of X-ray unobscured Type II AGN reported in 
Mainieri et al. (2007). 

An interesting source is XID=#2608 which has b een c lassi- 
fied as a Compton-T hick AGN by iMainieri etaT] d2007l) and 
lHasingere t al. ( 2007 ) but its optical s pectru m is that of an emis- 
sion line galaxy. In lHasinger et alj (|2007) a small number of 
sources (including XID=#2608) was found to have hardness ra- 
tios that could be interpreted as being due to heavily absorbed 
(possibly Compton thick) high energy spectra with some frac- 
tion of leaking unabsorbed soft flux. The solid green line in Fig. 
QT]represents the expected tracks occupied by leaking Compton 
thick sources in the HR1-HR2 planqj at z=0. The line has 
been computed with a pure reflection model with a fraction of 
1%,3%,10% and 30% (from top to bottom) of the flux from the 
central source leaking out. 

In particular the source ID=#2608 shows X-ray colours typi- 
cal of a spectrum dominated by a pure reflection component 
with ~3% of the original flux leaking out. Another source, 
XID=#131, shows X-ray colours consistent with Compton-thick 
AGN with a small fraction of leaking flux. We note that a 1% 
fraction of Compton-thick AGN is consistent with the predic- 
tions of XRB models at the flux limit of this subsample (i.e. 2- 



10 keV flux >10 erg cm s ) and with the source counts of 



9 The position in this HR1-HR2 plane of the tr ack of leaking 
Comp ton thick objects is different from that shown in Hasing er et al.1 
d2007l) because of the difference in the energy range of the hard band 
(i.e. 2-8 keV here, 2-4.5 keV in Hasin ger et all d2007h ). 
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Compton-thick obje cts measured in a collection of surveys by 
iBrunneret al1(l2008l) . 

The objects classified as emission and absorption line galax- 
ies are spread over the entire luminosity-hardness ratio plane 
(see Fig. IT~2T > and their nature can be explained as a mixture 
of star forming galaxies, Type II AGN and XBONGs (see e.g . 
Comastri et aljliooalCaccianiga et al.ll2007UCivano et al.ll2007t 
Cocchia et al.ll2007l) . A more detailed analysis of their multi- 
wavelength properties will be the subject of a forthcoming pub- 
lication. 



6. Summary 

In this paper we presented a pointlike source catalogue in the 
XMM-COSMOS survey. The survey covers an area of 2.13 deg 2 
in the equatorial sky. The field has been observed with 55 XMM- 
Newton pointings for a total exposure time of ~1.5 Ms. We 
achieved an almost uniform exposure of ~40 ks on the field. 
We detected a total number of 1 62 1 , 1111 and 25 1 sources in the 
0.5-2 keV, 2-10 keV and 5-10 keV energy band, respectively, 
for a total of 1887 independent sources detected with det_ml>10 
in at least one band. The survey has a limiting flux of ~1.7xl0" 15 
erg crrT 2 s _1 , ~9.3xl0~ 15 erg cirT 2 s and ~1.3xl0~ 14 erg 
cm" 2 s" 1 in the 0.5-2 keV, 2-10 keV and 5-10 keV energy 
band, over 90% of the area. 

Together with the source catalogue we derived logN-logS rela- 
tions with high statistics in the flux interval sampled by the sur- 
vey. The logN-logS relations are in good agreement with most of 
the X-ray surveys published in the literature. We compared our 
source c ounts with the most recent XRB popu lation synthesis 
models dGilli et al]l2007t iTreister & Urrvi r2006) and found that 
they agree within 10% with our data in the 5-10 keV and 2-10 
keV energy bands. In the 0.5-2 keV band both models deviate 
from the XMM-COSMOS data by about 10%-30% suggesting, 
that further improvements in the modeling are required. We iso- 
lated a subsample of X-ray bright sources for which optical spec- 
troscopy is available. About 65% of them have optical and X-ray 
properties typical of Type I AGN and ~ 15% of Type II AGN. In 
the subsample of sources with a good optical spectrum and good 
counting statistics, the number of candidate Compton thick (1-2) 
AGN is fully consistent with the expectations of XRB population 
synthesis models. By combining X-ray colours and optical spec- 
troscopy we found that 20% of the sources do not show, in the 
optical band, evident signatures of AGN activity although their 
X-ray luminosities are typical of AGN. Additonall y, we consider 
XMM -COSMOS as a pathfinder for the eROSITA dPredehl et alj 
2006) X-ray telescope which will be launched in 2012 and that 
will perform an all sky survey with sensitivities comparable to 
those presented here. 
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36.* 

41. n. 
43.61 

49.0 

58. as 

47. (i| 
40. S 
45.(3 

40. ©. 

41. * 
48.0^. 

58.H 
50.ll 
51.4 
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45 §i 

39.(& 
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53.34 
53.24 
42.46 



" : Estimated in the 2-8 keV band. 
b : Estimated in the 4.5-10 keV band. 
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Table 4. Source number counts. 



Log(S) 


Nt>S , l 




Nt>S , l 


Area 


Ni>S , l 




-i -i 
erg cm s 


deg 


deg z 


deg 


deg^ 


deg 


deg^ 




0.5-2 keV 




2-10 keV 




5-10 keV 




-12.8 


1.51 +0.81 


2.13 


3.86 ±1.33 


2.13 


1.04 ±0.66 


2.13 


-12.9 


2.45 ±1.05 


2.13 


6.20 ±1.69 


2.13 


1.98 ±0.94 


2.13 


-13.0 


3.39 ±1.24 


2.13 


8.55 ±1.99 


2.13 


2.92 ±1.15 


2.13 


-13.1 


5.26 ±1.56 


2.13 


11.84 ±2.35 


2.13 


4.79 ±1.48 


2.13 


-13.2 


8.08 ±1.94 


2.13 


13.25 ±2.48 


2.13 


7.14 ±1.82 


2.13 


-13.3 


10.43 ±2.20 


2.13 


24.51 ±3.39 


2.13 


8.55 ±1.99 


2.13 


-13.4 


14.65 ±2.61 


2.13 


32.96 ±3.93 


2.13 


11.84 ±2.35 


2.13 


-13.5 


19.35 ±3.01 


2.13 


50.33 ±4.86 


2.13 


17.94 ±2.89 


2.13 


-13.6 


25.92 ±3.48 


2.13 


77.56 ±6.03 


2.13 


24.98 ±3.42 


2.13 


-13.7 


40.48 ±4.35 


2.13 


108.07 ±7.12 


2.13 


39.07 ±4.28 


2.13 


-13.8 


56.44 ±5.14 


2.13 


150.33 ±8.40 


2.13 


55.57 ±5.11 


2.10 


-13.9 


76.63 ±5.99 


2.13 


208.52 ±9.90 


2.09 


84.72 ±6.36 


1.82 


-14.0 


102.45 ±6.93 


2.13 


277.20 ±11.48 


1.98 


123.75 ±8.26 


0.95 


-14.1 


131.09 ±7.84 


2.13 


361.50 ±13.31 


1.77 


166.26 ±13.33 


0.23 


-14.2 


166.77 ±8.85 


2.13 


491.38 ±16.27 


1.33 


212.88 ±43.57 


0.02 


-14.3 


217.00 ±10.09 


2.13 


620.89 ±19.91 


0.67 






-14.4 


273.34 ±11.33 


2.13 


766.57 ±29.63 


0.14 






-14.5 


324.98 ±12.35 


2.13 


984.00 ±106.76 


0.01 






-14.6 


398.23 ±13.67 


2.13 










-14.7 


480.54 ±15.04 


2.06 










-14.8 


581.69 ±16.68 


1.80 










-14.9 


713.44 ±19.05 


1.37 










-15.0 


842.39 ±21.69 


1.00 










-15.1 


930.12 ±24.30 


0.44 










-15.2 


1027.83 ±35.07 


0.08 










-15.3 


1201.41 ±177.09 


0.01 











